hep-ph/9911273 

CERN-TH/99-317 
RM3-TH/99-11 
Edinburgh 99/18 



Resummation of Singlet Parton Evolution at Small x 

Guido Altarelli,"' 6 Richard D. Ball '- Q and Stefano Forte c Q 

"Theory Division, CERN 
CH-1211 Geneva 23, Switzerland 

6 Universita di Roma Tre 
and 

C INFN, Sezione di Roma Tre 
Via della Vasca Navale 84, 1-00146 Rome, Italy 

^Department of Physics and Astronomy, University of Edinburgh, 
Mayfield Road, Edinburgh EH9 3JZ, Scotland 

Abstract 

We propose an improvement of the splitting functions at small x which overcomes the apparent 
problems encountered by the BFKL approach. We obtain a stable expansion for the x-evolution 
function x{M) near M = by including in it a sequence of terms derived from the one- and two- 
loop anomalous dimension 7. The requirement of momentum conservation is always satisfied. 
The residual ambiguity on the splitting functions is effectively parameterized in terms of the 
value of A, which fixes the small x asymptotic behaviour x~ x of the singlet parton distributions. 
We derive from this improved evolution function an expansion of the splitting function which 
leads to good apparent convergence, and to a description of scaling violations valid both at 
large and small x. 
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1. The theory of scaling violations in deep inelastic scattering is one of the most solid 
consequences of asymptotic freedom and provides a set of fundamental tests of QCD. At large 
Q 2 and not too small but fixed x the QCD evolution equations for parton densities provide 
the basic framework for the description of scaling violations. The complete splitting functions 
have been computed in perturbation theory at order a s (LO approximation) and a 2 (NLO) @. 
For the first few moments the anomalous dimensions at order are also known ||. 

At sufficiently small x the approximation of the splitting functions based on the first few 
terms in the expansion in powers of a s is not in general a good approximation. If not for 
other reasons, as soon as x is small enough that a s ^ ~ 1, with £ = logl/x, all terms of order 
a s (a s £) n and a 2 (a s £,) n which are present f§] in the splitting functions must be considered in 
order to achieve an accuracy up to terms of order af. In terms of the anomalous dimension 
j(N,a s ), defined as the iV-th Mellin moment of the singlet splitting function (actually the 
eigenvector with largest eigenvalue), these terms correspond to sequences of the form (a s /N) n 
or a s (a s /N) n . In most of the kinematic region of HERA || the condition a s £ ~ 1 is indeed true. 
Hence, in principle one could expect to see in the data indications of important corrections to the 
approximation || [7j of splitting functions computed only up to order a 2 and the corresponding 
small x behaviour. In reality this appears not to be the case: the data can be fitted quite 
well by the evolution equations in the NLO approximation |7], |j. Of course it may be that 
some corrections exist but they are hidden in a redefinition of the gluon, which is the dominant 
parton density at small x. While the data do not support the presence of large corrections in 
the HERA kinematic region || the evaluation of the higher order corrections at small x to the 
singlet splitting function from the BFKL theory [ |TU| , |TT| , |T2[ appears to fail. The results of the 
recent calculation |13[ |14], [HJ of the NLO term xi of the BFKL function x = c^sXa + ^Xi--- show 
that the expansion is very badly behaved, with the non leading term completely overthrowing 
the main features of the leading term. Taken at face value, these results appear to hint at 
very large corrections to the singlet splitting function at small x in the region explored by 
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In this article we address this problem and propose a procedure to construct a meaningful 
improvement of the singlet splitting function at small x, using the information from the BFKL 
function. We start by defining an alternative expansion for the BFKL function x(M) which, 
unlike the usual expansion, is well behaved and stable when going from LO to NLO, at least 
for values far from M — 1. This is obtained by adding suitable sequences of terms of the 
form (a s /M) n or a s (a s /M) n to a s xo or a 2 xi respectively. The coefficients are determined by 
the known form of the singlet anomalous dimension at one and two loops. This amounts to a 
resummation |17| of (a s log Q 2 / fi 2 ) n terms in the inverse M-Mellin transform space. This way of 
improving x is completely analogous to the usual way of improving 7 |TJ|. One important point, 
which is naturally reproduced with good accuracy by the above procedure, is the observation 
that the value of x(M) at M = is fixed by momentum conservation to be x(0) = 1. This 
observation plays a crucial role in formulating the novel expansion and explains why the normal 
BFKL expansion is so unstable near M = 0, with xo ~ ^/M, xi ~ —1/M 2 and so on. This 
rather model-independent step is already sufficient to show that no catastrophic deviations 
from the NLO approximation of the evolution equations are to be expected. The next step is 
to use this novel expansion of x t° determine small x resummation corrections to add to the LO 
and NLO anomalous dimensions 7. Defining A as the minimum value of x, x(.M min ) = A, and 
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using the results of ref. [19], a meaningful expansion for the improved anomalous dimension 
is written down in terms of Xo, Xii an d A. The large negative correction to A /a s = Xo(l/2) 
induced by a s Xi, that is formally of order a s but actually is of order one for the relevant values 
of a s , suggests that A should be reinterpreted as a nonperturbative parameter. We conclude 
by showing that the very good agreement of the data with the NLO evolution equation can be 
obtained by choosing a small value of A, compatible with zero. 

2. We consider the singlet parton density 

G(£, t) = x[g(x, Q 2 ) + k q ® q(x, Q% (1) 

where £ = logl/x, t = logQ 2 //z 2 , g(x,Q 2 ) and q(x,Q 2 ) are the gluon and singlet quark parton 
densities, respectively, and k q is such that, for each moment 

poo 

G(N,t) = / d{e- N tG({,t), (2) 
Jo 

the associated anomalous dimension j(N,a s (t)) corresponds to the largest eigenvalue in the 
singlet sector. At large t and fixed £ the evolution equation in iV-moment space is then 

jG(N,t)= 1 (N,a s (t))G(N,t), (3) 

where a s (t) is the running coupling. The anomalous dimension is completely known at one and 
two loop level: 

7 (iV, a,) = a s j {N) + a 2 s ^(N) + ... . (4) 

As j(N, a s ) is, for each N, the largest eigenvalue in the singlet sector, momentum conservation 
order by order in a s implies that 

7 (l,a.)=7o(l) = 7i(l) = - = 0. (5) 
Similarly, at large £ and fixed Q 2 , the following evolution equation for M moments is valid: 

^G(£,M)= X (M,a a )G(£,M), (6) 

where _ 

POO 

G(i,M)= / dte- Mt G(U), (7) 
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and xiM, a s ) is the BFKL function which is now known at NLO accuracy: Q 

X (M,a s ) = a sXo (M) + a 2 sX i(M) + . . . . (8) 

In eq. (g) the coupling a s is fixed. The inclusion of running effects in the BFKL theory 
is a delicate point. To next-to-leading order in a s (i.e. to NLLx), running effects can be 
included |L9|, |2(J by adding to Xi a term proportional to the first coefficient /3 = yn c — |n/ 

^^Note that the normalization conventions for xo ancl Xl used here are different from those used in either of 
refs.ra,|T9|. 
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Figure 1: Graphical representation of different expansions of (a) 7 and (b) x i n powers of a s and 
1/N (a) and 1/M (b) to order m and n respectively, and the different relations between these 
expansions. Vertical lines correspond to terms of the same fixed order in a s : for example the one 
loop anomalous dimension 70 contains terms with m = 1, n = 1,0,— 1,— 2,.... Diagonal lines 
correspond to terms of the same order in a s at fixed ^ (a) or || (b): for example •y s (a s /N) 

contains terms with m — n — 1, 2, 3, The sum of terms in a vertical line of the 7 plot 

is related by duality eqs. (|10|, lTD to the sum of terms in a diagonal line in the x plot an d 
conversely (marked by the same line style). The solid lines denote terms of the same order in 
the "envelope" or "double leading" expansion discussed in the text. 



of the /3-function. Since furthermore the extra term depends on the definition of the gluon 
density, it is also necessary to specify the choice of factorization scheme: here we choose the 
MS scheme, so that the xi that we will consider in the sequel is given by [|l]| 

Xi (M) = ^n 2 J(M) + ^p n c ((2^'(l) - ip'(M) - ^(1 - M)) + ^Xo(M) 2 , (9) 



where the function 5 is defined in the first of ref. 13 



In the region where Q 2 and 1/x are both large the t and £ evolution equations, i.e. eqs.(||,|). 
are simultaneously valid, and their mutual consistency requires the validity of the "duality" 
relation 0, 0]: 

X (l(N,a s ),a s )=N, (10) 

and its inverse 

1 ( X (M,a s ),a s )=M. (11) 

Using eq. (|T0D, knowledge of the expansion eq. (§) of x(M,a s ) to LO and NLO in a s at fixed 
M determines the coefficients of the expansion of 7 (AT, a s ) in powers of a s at fixed jf: 

7 (JV, a s ) = 7s (f ) + a slss (f) +■■■, (12) 
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where 7 S and 7 SS contain respectively sums of all the leading and subleading singularities of 7 
(see fig. 1), 

Xo(7 s (t») = (13) 

This corresponds to an expansion of the splitting function in logarithms of x: if for example 
we write 

= £/(£)' (15) 

(where g± = n c /ir, g^ = #3 = ,g± = 2£(3)(n c /7r) 4 , . . .), then the associated splitting 
function 

r c+ioo dN 

and similarly for the subleading singularities P ss (q s {), etc. 

Likewise, the inverse duality eq. (PD relates the fixed order expansion eq. (§]) of j(N, a s ) 
to an expansion of x(M, a s ) in powers of a s with ^ fixed: if 

X (M, a,) = x.{%) + a sXss (^) + ■■-, (17) 

where now Xs(ff) an d Xss(ff-) contain the leading and subleading singularities respectively of 
x(M,a s ), then 

7o(x*(§ )) = -, (18) 

Ot s 

Xss("m) = ~~ ~> T , Ta, \\ ■ (1^) 



rc+ioo AN 00 

p ^ b U ^^^^ra^' (16) 



7ife(f)) 
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In principle, since xo and xi are known, they can be used to construct an improvement of 
the splitting function which includes a summation of leading and subleading logarithms of x. 



However, as is now well known, the calculation [13, 14, 15| of Xi has shown that this procedure 



is confronted with serious problems. The fixed order expansion eq. (§) is very badly behaved: 
at relevant values of a s the NLO term completely overwhelms the LO term. In particular, near 
M — 0, the behaviour is unstable, with xo ~ xi ~ — 1/M 2 . Also, the value of x near 

the minimum is subject to a large negative NLO correction, which turns the minimum into a 
maximum and can even reverse the sign of x a t the minimum. Finally, if one considers the 
resulting 7 S and 7 SS or their Mellin transforms P s (x) and P ss (x) one finds that the NLO terms 
become much larger than the LO terms and negative in the region of relevance for the HERA 
data [TBJ. We now discuss our proposals to deal with all these problems. 



Our first observation is that a much more stable expansion for xiM) can De obtained if we 
make appropriate use of the additional information which is contained in the one and two loop 
anomalous dimensions 70 and 71. Instead of trying to improve the fixed order expansion eq. (KJ) 
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of 7 by all order summation of singularities deduced from the fixed order expansion eq. (§) of 
X, we attempt the converse: we improve Xa(M) by adding to it the all order summation of 
singularities Xs eq. (|18D deduced from 70, Xi(^0 by adding to it Xss deduced from ji eq. (|19"D, 
and so on. It can then be seen that the instability at M = of the usual fixed order expansion of 
X was inevitable: momentum conservation for the anomalous dimension, eq. (JD, implies, given 
the duality relation, that the value of x(M) at M = is fixed to unity, since from eq. (10) we 
see that at N = 1 

X(7(l,a.),a.)=x(0,a a ) = l. (20) 

It follows that the fixed order expansion of x must be poorly behaved near M = 0: a simple 
model of this behaviour is to think of replacing a s /M with a s / (M + a s ) = a s /M — a 2 /M 2 + ... 
in order to satisfy the momentum conservation constraint. We thus propose a reorganization of 
the expansion of x a "double leading" (DL) expansion, organized in terms of "envelopes" 
of the contributions summarized in fig. lb: each order contains a "vertical" sequence of terms 
of fixed order in a s , supplemented by a "diagonal" resummation of singular terms of the same 
order in a s if a s /M is considered fixed. To NLO the new expansion is thus 



X (M,a s ) = [a sX o(M)+ Xs (§) 



■kM 



+a s [a sX x{M) + Xss (f ) ~ <*.(■& + £)- fo] + - (21) 

where the LO and NLO terms are contained in the respective square brackets. Thus the LO 
term contains three contributions: Xo(M) is the leading BFKL function eq. (§), Xs{oi s /M) 
eq. fll8 ) are resummed leading singularities deduced from the one loop anomalous dimension, 
and n c a s / '(jrM) is subtracted to avoid double counting. At LO the momentum conservation 
constraint eq. (p0|) is satisfied exactly because 70 (1) = and [xo(M) — ~ M 2 near M — 0. 
At NLO there are again three types of contributions: Xi(M) from the NLO fixed order calcu- 
lation (eq. (H)), the resummed sub leading singularities Xss(o^ s /M) deduced from the two loop 
anomalous dimension, and three double counting terms, fo — 0, fi — — n/(13+ 10?i|)/(367r 2 r^) 
and / 2 = nl(ll+2n f / nl) / (12n 2 ) (corresponding to those terms with (m, n) = (1, 0), (2, 1), (2, 2) 
respectively in fig. lb). Note that at the next-to-leading level the momentum conservation con- 
straint is not exactly satisfied because the constant contribution to Xi does not vanish in MS, 
even though it is numerically very small (see fig. 2). It could be made exactly zero by a re- 
finement of the double counting subtraction but we leave further discussion of this point for 
later. 

Plots of the various LO and NLO approximations to x are shown in fig. 2. In this and 
other plots in this paper we take a s = 0.2, which is a typical value in the HERA region, and 
the number of active flavours rif = 4. We see that, as discussed above, the usual fixed order 
expansion eq. in terms of Xo an d Xi is very unstable. However, the new expansion eq. (|21"D 
is stable up to M ^ 0.3 — 0.4. Furthermore, in this region, x evaluated in the double leading 
expansion (^) is very close to the resummations of leading and subleading singularities eq. (|PTD 
obtained by duality eq. ( |Il||l9"l) from the one and two loop anomalous dimensions. This shows 



that in this region the dominant contribution to x> an d thus to 7 , comes from the resummation 
of logarithms of Q 2 /fi 2 with Q 2 > fi 2 . 

Beyond M ~ 0.4, the size of the contributions from collinear singular and nonsingular terms 
becomes comparable (after all here Q 2 ~ /z 2 ), but the calculation of the latter (from the fixed 
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Figure 2: Plots of different approximations to the x function discussed in the text: the BFKL 
leading and next-to- leading order functions eq. (||), a s xo and a s xo + o? s Xi (dashed); the LO and 
NLO dual a s Xs and a s x s + Oi 2 s x ss of the one and two loop anomalous dimensions (solid), and 
the double- leading functions at LO and NLO defined in eq. ( pTI ) (dotdashed). Note that the 
double-leading curves coincide with the resummed ones at small M, and with the fixed order 
ones at large M. 



expansion eq. (§)) has become unstable due to the influence of the singularities at M = 1. 
No complete and reliable description of x seems possible without some sort of stabilization of 
these singularities. However, since they correspond to infrared singularities of the BFKL kernel 
(specifically logarithms of Q 2 /[i 2 with Q 2 <C /i 2 ) this would necessarily be model dependent. In 
particular, such a stabilization cannot be easily deduced from the resummation of the M = 
singularities: the original symmetry of the gluon-gluon amplitude at large s is spoiled by 
running coupling effects and by unknown effects from the coefficient function through which it 
is related to the deep-inelastic structure functions, in a way which is very difficult to control 
near the photoproduction limit M — > 1. We thus prefer not to enter into this problem: rather 
we will discuss later a practical procedure to bypass it. 

The results summarized in fig. 2 clearly illustrate the superiority of the new double leading 
expansion of x over the fixed order expansion, and already indicate that the complete x function 
could after all lead to only small departures from ordinary two loop evolution. 

3. Having constructed a more satisfactory expansion eq. (|2l|) of the kernel Xi we now derive 
from it an improved form of the anomalous dimension 7 to be used in the evolution eq. (|3|), in 
order to achieve a more complete description of scaling violations valid both at large and small 
x. In principle, this can be done by using the duality relation eq. (|TUD, which simply gives the 
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function 7 as the inverse of the function x- However, in order to derive an analytic expression 
for 7(iV, a s ) which also allows us to clarify the relation to previous attempts we start from the 
naive double-leading expansion of 7 |L8| in which terms are organized into "envelopes" of the 
contributions summarized in fig. la in an analogous way to the double leading expansion ( |2T| ) 
of x- 



N 



a*7o(A0+7s(f) - 

a s n(N)+ lss (f) 



ry ( M, -L 21 
U S I TV 2 



N2 + N I 



(22) 



where now e 2 = g£ = 0, e 1 = gi = n f n c (5 + 13 / (2H 2 )) / (18ir 2 ) and e = -(^nl + n^/iQim 2 ). 
In this equation, the leading and subleading singularities 7 S and 7 SS are obtained using duality 
eq. ([H]) from xo and Xii an d summed up to give expressions which are exact at NLLx. These 
are then added to the usual one and two loop contributions, and the subtractions take care of 
the double counting of singular terms. 

It can be shown that the dual of the double leading expansion of x ec i- (HI) coincides with 
this double leading expansion of 7 eq. (|22|) order by order in perturbation theory, up to terms 
which are higher order in the sense of the double leading expansions. However, it is clear that 
these additional subleading terms must be numerically important. Indeed, it is well know that 
at small N the anomalous dimension in the small- a; expansion eq. (|T2|) is completely dominated 
by lss{oi s /N) which grows very large and negative, leading to completely unphysical results 
in the HERA region |16[|. It is clear that this perturbative instability will also be a problem 



in the double leading expansion eq. (^). On the other hand, we know from fig. 2 that the 
exact dual of x m double leading expansion is stable, and not too far from the usual two loop 
result. The origin of this instability problem, and a suitable reorganization of the perturbative 
expansion which allows the resummation of the dominant part of the subleading terms have 
been discussed in ref. [ 19| . After this resummation, the resulting expression for 7 in double 



leading expansion will be very close to the exact dual of the corresponding expansion of x- 



The procedure of ref. |19| can be interpreted in a simple way whenever the all-order "true" 
function x^M-, a s) possesses a minimum at a real value of M, M min , with < M min < 1 
(although the final result for the anomalous dimension will retain its validity even in the absence 
of such minimum). Using A to denote this minimum value of Xi 



A = X (M„ 



a. 



A + AA, 



Ao = a s Xo 



4n c 



a q In 2. 



(23) 



The instability turns out to be due to the fact that higher order contributions to 7 must 
change the asymptotic small x behaviour from x~ x ° to x~ x . The starting point of the proposed 
procedure consists of absorbing the value of the correction to the value of x a ^ the minimum 
into the leading order term in the expansion of x'- 



X (M,a s ) = a sXo (M) + a 2 sX i(M) + . . . 

= (a sX o(M) + AA) + a 2 s xi(M) + ■ 



(24) 



where x n (M) = Xn(M) — c n , with c n chosen so that x n (M) no longer leads to an 0(a") shift 
in the minimum. Since the position M m « n of the all-order minimum is not known, one must 
in practice expand it in powers of a s around the leading order value M = ^, so at higher 
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orders the expressions for the subtraction constants c n can become quite complicated functions 



of Xi an d their derivatives at M — | |T9]. However at NLO we have simply c\ = Xi(f)j so 
AA = a 2 sX x{\) + • • •■ 

A stable expansion of 7 in resummed leading and subleading singularities can now be ob- 
tained from the duality eqs. (p~3| , p^4| ,. . . ) by treating xo + AA as the LO contribution to Xi an d 
the subsequent terms x% as perturbative corrections to it. Of course, since the reorganization 
eq. fl24|) amounts to a reshuffling of perturbative orders, to any finite order the anomalous 
dimension obtained in this way will be equal to the old one up to formally subleading correc- 
tions. Explicitly, we find in place of the previous expansion in sums of singularities eq. (|T2| ) the 
resummed expansion 

7 (JV, a,) = is (^u) + a s % s (j^) +..., (25) 

where 



7ss (jv_aa) — ^ ss (jv-aa) + , ( 

Xo(7- 



,JV-AA 

In terms of splitting functions this resummed expansion is simply 



(26) 



xP(x,a s ) = a s e ?AA [P s (a s O + a s P ss (a s Q + . . .] 

= a s e^ AX [P s (a s O + a s P ss (a s O ~ ^Psi^sO + ■■■}■ (27) 



The expansion is now stable in the sense that it may be shown that P ss (a s ^) / P s (a s £) 
remains bounded as £ ^00: subleading corrections will then be small provided only that a s 
is sufficiently small. This result may be shown to be true to all orders in perturbation theory, 
using an inductive argument. 

We can thus replace the unresummed singularities 7 S and 7 SS in eq. (|2"2"D with the resummed 
singularities eq. (f25"|) to obtain a double leading expansion with stable small x behaviour: 



7 (iV,a s ) = |a s7o (A0 + 7s (^^) - a s ^ 



a s7l (N) + lUj^fkx) ~ a s (jk + f ) - e + • • • . (28 



Momentum conservation is violated by the resummation because 7^ and 7 SS and the subtraction 
terms do not vanish at N = 1. It can be restored by simply adding to the constant eo a further 
series of constant terms beginning at 0{a 2 s ): these are all formally subleading in the double 
leading expansion. This constant shift in 7 is precisely analogous to the shift made on x hi 
eq. (p3|) which generated the resummation. 

It is important to recognize that there is inevitably an ambiguity in the double counting 
subtraction terms in eq. (|28|). For example, at the leading order of the double leading expansion 
instead of subtracting ^f- we could have subtracted -^jfr&xj, since this differs only by formally 
subleading terms: A A = 0(a 2 ), so 
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Following the same type of subtraction at NLO, the resummed double leading anomalous 
dimension may thus be written as 



a sll (N) + lss(A- x ) + - + A) - eoj + ■ • • (30) 

3 

The extra term at NLO comes from the first correction in eq. (ffif) , which is of order and 



N 2 , 

thus a subleading singularity. The characteristic feature of this alternative resummation is 
that the fixed order anomalous dimensions 70, 71 are preserved in their entirety, including the 
position of their singularities. As with the previous expansion eq. momentum conservation 
may be imposed by adding to eo a series of terms constant in N and starting at 0(a 2 s ). 

This completes our procedure of inclusion of the most important part of the subleading 
corrections, as we shall see shortly by a direct comparison of the resummed expansions eq. fl28|) 
and eq. ([3~0"1) with the exact dual of x evaluated according to eq. (|2l|) . In the sequel we will 
discuss the phenomenology based on the two resummed expansions eq. (^8|) and eq. fl30"D on 
an equal footing, taking the spread of the results as an indication of the residual ambiguity 
due to subleading terms. Although formally the differences between the two expansions are 
subleading, we will find that in practice they may be quite substantial, because AA may be 
large. 

4. So far we have constructed resummations of the anomalous dimension and splitting function 
which satisfy the elementary requirements of perturbative stability and momentum conserva- 
tion. This construction relies necessarily on the value A of x near its minimum, since it is this 
which determines the small x behaviour of successive approximations to the splitting function. 
In order to obtain a formulation that can be of practical use for actual phenomenology, we will 
need however to improve the description of x{M) in the "central region" near its minimum 
M m i n , since as we already observed, we cannot reliably determine the position and value of 
the minimum of \ without a stabilization of the M = 1 singularity. Indeed, we can see from 
fig. 2 that in the central region \ evaluated in the double leading expansion is dominated by 
the presumably unphysical M = 1 poles of Xi an d at NLO this means that it actually has no 
minimum, becoming rapidly negative. However, one can use the value A of the true x at the 
minimum as a useful parameter for an effective description of the x function around M = 1/2. 
Indeed, AA as estimated from its next-to-leading order value a^Xi(l/2) turns out to be of the 
same order as A for plausible values of a s , a feature which can be also directly seen from 
fig. 2. This supports the idea that A and AA are not truly perturbative quantities: in general 
we expect that the overall shift of the minimum will still be of the order of Aq and negative. 



It is this order transmutation that makes the impact of the resummations eq. (|28| , |30"1) , and the 
differences between them, quite substantial. 

In fig. 3 and fig. 4 we display the results for the resummed anomalous dimensions in the two 
different expansions eq. (|^) and eq. (0) respectively, each computed at next-to-leading order. 
In both figures we show for comparison the fixed order anomalous dimension a s 7o(A^) + a^7i(A^) 
eq. ([|). Also for comparison, we show the exact dual of x computed at NLO in the double 
leading expansion eq. obtained from eq. ( [TD|) by exact numerical inversion. This curve is 
thus simply the inverse of the corresponding curve already shown in fig. 2. 
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Figure 3: Comparison of the anomalous dimension 7 evaluated at NLO in the resummed 
expansion eq. (^3) for three different values of A (dashed) with the usual fixed order perturbative 
anomalous dimension (also at NLO) eq. (f|) (dotted) and that obtained by exact duality from 
X at NLO in the expansion eq. (^T|) as displayed in fig. 2 (solid). The unresummed 7 eq. (|22|) 
is also shown at NLO. Notice that the A = 0.21 curve is very close to the two loop anomalous 
dimension down to the branch point at iV = A. 



In fig. 3 we show the anomalous dimension computed at NLO using the resummation 
eq. (p8|), for A = Ao and A = 0. The first value corresponds to the LO approximation to A, while 
the second value is close to the NLO approximation when a s is in the region a s ~ 0.1 —0.2. We 
might expect the value of A as determined by the actual all-order minimum of x to lie within 
this range. Note that, in general, the resummed anomalous dimension has a cut starting at 
N = A, which corresponds to the x~ x power rise; for this reason our plots stop at this value of 
N. The A = curve, corresponding to the next-to-leading order approximation to A, is seen to 
be very close to the exact dual of x at NLO in the expansion eq. (^3), as already anticipated. 
This is to be contrasted with the corresponding unresummed anomalous dimension eq. (p2|) , 
which is also displayed in fig. 3, and is characterized by the rapid fall at small N discussed 
already in ref. |T6| , |19|j . This comparison demonstrates that indeed the perturbative reorgani- 
zation eliminates this pathological steep decrease. The resummed curve with A = and the 
exact dual of x become rather different for small N & 0.2. However, this is precisely the range 
of N which corresponds to the central region of M where we cannot trust the next-to-leading 
order determination of x- Finally, we show that we can choose a value of A ~ 0.21 such that 
the resummed anomalous dimension closely reproduces the two loop result down to the branch 
point at N = A. This shows that the absence of visible deviations from the usual two loop 
evolution can be accommodated by the resummed anomalous dimension. However this is not 



10 




0.5 1 1.5 2 2.5 



N 

Figure 4: Comparison of the anomalous dimension 7 evaluated at NLO in the resummed 
expansion eq. ( |5UD for two different values of A (dashed) with the usual fixed order perturbative 
anomalous dimension (also at NLO) eq. (f|) (solid) and that obtained by exact duality from \ 
at NLO as in fig. 3 (dotdash). Notice that the A = curve is virtually indistinguishable from 
the fixed order anomalous dimension for all values of N. 



necessarily the best option phenomenologically: perhaps the data could be better fitted by a 
different value of A if a suitable modification of the input parton distributions is introduced. 
It is nevertheless clear that large values of A such as A ~ Ao can be easily excluded within the 
framework of this resummation, since they would lead to sizeable deviations from the standard 
two loop scaling violations in the medium and large x region. 

The splitting functions corresponding to the anomalous dimensions of fig. 3 are displayed 
in fig. 5. The basic qualitative features are of course preserved: in particular, the curves with 
small values of A = and A = 0.21 are closest to the two loop result. However, on a more 
quantitative level, it is clear that anomalous dimensions which coincide in a certain range of N, 
but differ in other regions (such as very small N) may lead to splitting functions which differ 
over a considerable region in x. In particular, the A = 0.21 curve displays the predicted x~ x 
growth at sufficiently large £ (x ~ 10~ 4 ). The dip seen in the figure for intermediate values of £ 
is necessary in order to compensate this growth in such a way that the moments for moderate 
values of N remain unchanged. Note that the x~ x behaviour of the splitting function at small x 
is corrected by logs [pll: P s ~ £ -3 / 2 x~ A . If A = this logarithmic drop provides the dominant 

£ — ► OG 

large £ behaviour which appears in the figure. 

If the anomalous dimensions are instead resummed as in eq. (|50|), the results are as shown 
in fig. 4, again for the two very different values of A, A = and A = Ao- When A = the 
resummed anomalous dimension is now essentially indistinguishable from the two loop result. 
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Figure 5: The splitting functions corresponding to the anomalous dimensions of fig. 3. 



This is due to the fact that the simple poles at N = which are now retained in 70 and 71 
provide the dominant small N behaviour. The branch point at N = A in 7 S and 7 SS is then 
relatively subdominant. This remains of course true for all A < 0, and in practice also for small 
values of A such as A & 0.1. When instead A = Ao the result does not differ appreciably from 
the resummed anomalous dimension shown in fig. 3, since now the dominant small N behaviour 
is given by the branch point at iV = A , which is not affected by changes in the double counting 
prescription. Summarizing, the peculiar feature of the resummation eq. (|30D is that it leads to 
results which are extremely close to usual two loops for any value of A < 0, without the need 
for a fine-tuning of A. 

Finally, in fig. 6 we display the splitting functions obtained from the resummed anomalous 
dimensions of fig. 4. The A = Ao CcLSG IS, cLS expected, very close to the corresponding curve 
in fig. 5. However the A = curve is now in significantly better agreement with the two loop 
result than any of the resummed splitting functions of fig. 5, even that computed with the 
optimized value A = 0.21. Moreover, this agreement now holds in the entire range of £. This 
is due to the fact that the corresponding anomalous dimension is now very close to the fixed 
order one for all N > 0, and not only for N > A = 0.21. This demonstrates explicitly that one 
cannot exclude the possibility that the known small x corrections to splitting functions resum 
to a result which is essentially indistinguishable from the two-loop one. This however is only 
possible if A ~ 0. 

To summarise, we find that the known success of perturbative evolution, and in particular 
double asymptotic scaling at HERA can be accommodated within two distinct possibilities, 
both of which are compatible with our current knowledge of anomalous dimensions at small x, 
and in particular with the inclusion of corrections derived from the BFKL equation to usual 
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Figure 6: The splitting functions corresponding to the anomalous dimensions of fig. 4. 



perturbative evolution. One possibility, embodied by the resummed anomalous dimension 
eq. ( j30|) with A ^ 0, is that double scaling remains a very good approximation to perturbative 
evolution even if the x — > limit is taken at finite Q 2 . The other option, corresponding to the 
resummation eq. (^) with a small value of A, is that double scaling is a good approximation 
in a wide region at small x, including the HERA region, but eventually substantial deviations 
from it will show up at sufficiently small x. In the latter case, the best-fit parton distributions 
might be significantly differ from those determined at two loops even at the edge of the HERA 
kinematic region. Both resummations are however fully compatible with a smooth matching 



to Regge theory in the low Q region |22 



5. In conclusion, we have presented a procedure for the systematic improvement of the split- 
ting functions at small x which overcomes the difficulties of a straightforward implementation 
of the BFKL approach. The basic ingredients of our approach are the following. First, we 
achieve a stabilization of the perturbative expansion of the function x near M = through the 
resummation of all the LO and NLO collinear singularities derived from the known one- and 
two-loop anomalous dimensions. The resulting x function is regular at M — 0, and in fact, 
to a good accuracy, satisfies the requirement imposed by momentum conservation via duality. 
Then, we acknowledge that without a similar stabilization of the M = 1 singularity it is not 
possible to obtain an reliable determination of x i n the central region M ~ 1/2. However, we 
do not have an equally model-independent prescription to achieve this stabilization at M — 1. 
Nevertheless, the behaviour of x i n the central region can be effectively parameterized in terms 
of a single parameter A which fixes the asymptotic small x behaviour of the singlet parton 
distribution. This enables us to arrive at an analytic expression for the improved splitting 
function, which is valid both at small and large x and is free of perturbative instabilities. 
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This formulation can be directly confronted with the data, which ultimately will provide a 
determination of A along with a s and the input parton densities. The well known agreement of 
the small x data with the usual Q 2 evolution equations suggests that the optimal value of A will 
turn out to be small, and possibly even negative for the relevant value of a s . Such a value of 
A is theoretically attractive, because it corresponds to a structure function whose leading-twist 
component does not grow as a power of x in the Regge limit: it would thus be compatible 
with unitarity constraints, and with an extension of the region of applicability of perturbation 
theory towards this limit. 

Several alternative approaches to deal with the same problem through the resummation 
of various classes of formally subleading contributions have been recently presented in the 
literature. Specific proposals are based on making a particular choice of the renormalization 
scale |23) , or on a different identification of the large logs which are resummed by the £ evolution 
equation (|B|), either by a function of £ itself Plfl , or by a function of Q 2 |T7|, or both 



The main shortcoming of these approaches is their model dependence. For instance, in ref. |25] 
the value of A is calculated, and x is supposedly determined for all < M < 1. This however 
requires the introduction of a symmetrization of Xi which we consider to be strongly model 
dependent: indeed, in ref. [§5] it is recognized that their value of A only signals the limit of 
applicability of their computation. We contrast this situation with the approach to resummation 
presented here, which makes maximal use of all the available mo del- independent information, 
with a realistic parameterization of the remaining uncertainties. We expect further progress to 
be possible only on the basis of either genuinely nonperturbative input, or through a substantial 
extension of the standard perturbative domain. 
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